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ATTENUATION AND DISPERSION OF LONGITUDINAL WAVES 
IN VISCOUS, PARTLY IONIZED GASES* 

By David Kahn 
Electronics Research Center 


The propagation of longitudinal waves in three fluid, partly ionized gases is studied as a function 
of the collision and plasma-to- wave-frequency ratios. A collisional approach is taken with the full 
Boltzmann equations (together with Maxwell’s equations) serving as the starting point, the intention 
being to describe wave dispersion and absorption in the collisional regime, and to include the viscous 
and heat conduction effects on the wave propagation from the beginning as an integral part of the theory. 
While wave propagation in both the collisionless and continuum regimes has been extensively studied, 
collisional kinetic investigations have been more limited [1 J- [4 J and for the three fluid gas, the full 
Boltzmann equation approach employed in this paper for studying wave propagation has the advantage 
that all three fluid collision frequencies are used rather than only the neutral one of [ 3 ], and are obtained 
directly from the collision integrals. Since viscosity is a molecular interaction phenomenon, this direct 
identification with the collision integrals is consistent and desirable, the self collision frequencies de- 
termining the equilibration of the individual fluids and the interparticle collision frequencies determining 
the effective coupling between the fluids. On the other hand, the complicated equations that result from 
this approach lead to possible disadvantages: In order to simplify the equations we assume Maxwell 
molecule interactions involving the neutral particles and assume a long- and short-range Coulomb force 
(cut off at the Debye length) for interactions among the charged particles. While these are certainly ap- 
proximations to the true law of interaction, the error introduced is probably not very large as has been 
discussed previously [4], The equations are further simplified by assuming that the fluids are initially 
(prior to being perturbed) at the same temperature. This common temperature assumption leads to large 
ion wave damping because of phase mixing as the ion wave and particle speeds become comparable [lj. 
However, since this occurs predominantly in the Landau regime, for the purposes of this paper which is 
intended to determine the amount of collisional damping and collisional dispersion, the restriction ap- 
pears not to be too severe. 

For our collisional study we expand each of the velocity distribution functions, satisfying the cor- 
responding Boltzmann equation, in Sonine polynomials and solid harmonics (expansions about equilib- 
rium), and obtain the coupled set of equations for longitudinal waves shown below. The notations used 
are the 6 - functions <5 which are zero unless the interaction is between charged species in which case 
it is equal to unity and 6 g q = 1 if s refers to an electron or ion; otherwise it is zero; a propagation 
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constant k Q = KC nc /u normalized with respect to the adiabatic s-fluid sound speed C os -\| kg T 0 /m s ; 

5 ub 9 M\ U 6 

complex frequencies j n = (to* • n +icJl/ e i n^ where u n =0 and u e . are the usual electron and 
ion plasma frequencies and where v , and v ^ are the total electron, ion and neutral effect- 

ive cross-collision frequencies for momentum transfer, respectively; the effective self collision fre- 
quencies are included in the occurring in the higher order equations (3) and (4) below. Finally 

p = m 0 n n is the equilibrium mass density of the s-fluid, s denoting the electron, ion, or neutral 
os s os 

particles: 
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From this set of 12 equations, 3 each for the momentum (Aq]), energy (Ajq), stress (Aq 2 ) and heat 
flow (A^i), we immediately form the longitudinal wave dispersion relation as the determinantal equation 
DET A = 0 with the elements of this determinant being just the coefficients of the Ap^ g terms. The dis- 
persion relation has been solved numerically for the propagation constant over a large range of collision, 
plasma, and wave frequency ratios for several different plasmas ranging from slightly ionized to fully 
ionized, and including several plasmas which are intended to approximate some experimental situations 
1 5] . Because of limited space, we show here only the results of a couple of cases obtained: The neutral 
wave in a slightly ionized gas and the ion wave in a fully ionized gas (figures 1-4). These examples, 
though picked for their relative simplicity, do, however, illustrate the effects of collisional coupling on 
the wave propagation and the general result of a maximum in the absorption constant when the self col- 
lision frequency is approximately equal to the wave frequency. Further information is obtained by noting 
several limiting cases that were taken. In the limit of very large self-collision-to-wave-frequency ratios, 
the dispersion relation reduces to the inviscid dispersion relations of {4] and, in the fully ionized case, 
of [ 6 J , since the viscosity vanishes in this limit where i/ ss /w — - oo. If we also idealize to the case of 
negligible interparticle collision frequencies in this inviscid limit, we obtain just the usual nondissipa- 
tjive adiabatic dispersion relations discussed in (2] . 

In all cases solution of the dispersion relation yields the six roots C 0S /C I = 0.606 and 1.59 when 
the collision-to-wave-frequency ratios are sufficiently small (a result also obtained for neutral gases by 




Figure 2. - Neutral particle wave absorption 
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Grad [ 7] in this limit); these roots will serve as convenient labels for the waves picked out for discus- 
sion here. In figures 1 and 2 we show the C^ n /C = 0.606 root for a slightly ionized gas (the neutral to 

ionized particle density was taken as 2x10^) as a function of the neutral- neutral collision-to-wave- 
frequency ratio v nn /u . We find an absorption peak around tv l where the individual collisions 

are most effective in extracting energy from the wave; as v nn /u) (and also v j n / u ) further in creases to 
values greater than unity, the wave speed approaches the total adiabatic sound speed C ot = ^| /P 0 t 
appropriate for the slightly ionized collision-dominated gas. 

In figures 3 and 4 we show both ion wave roots C Q j /C = 0.606 and 1.59 of the fully ionized ex- 
ample plotted as a function of the ion-electron collision to wave frequency ratio ^- e /u (or, as indicated 
in the figures, as a function of and Wj /•*> , as well). The roots exhibit marked differences in 

their dependence on these ratios. The slow root (1.59) has a strong phase speed dependence on w j 
and becomes strongly damped just as the collisional regime becomes established (v.. /u> > l). The fast 
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root, on the other hand, shows no phase speed dependence upon the ion plasma frequency w j, maintaining 
a constant speed until v ^ /u « 1 at which point the wave absorption shows its characteristic maximum 
and the phase speed decreases to a value near the adiabatic ion sound speed. As the ions and electrons 
become coupled through collisions for i/j e / u > 1. the absorption correctly approaches zero and the 
wave travels with the adiabatic sound speed as it should in the collision-coupled gas. 
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